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We consider the factorization of hadronic physical observables with transverse momentum depen- 
dence when two collinear sectors are involved and which interact through emission of soft gluons. 
By proving a novel property of the relevant soft function as a geometric mean of two "hemisphere" 
soft functions, one can combine those functions with the two collinear ones to obtain a modified fac- 
torization theorem in which all quantities are well-defined and each "generalized" collinear quantity 
belongs to distinct region in rapidity space. Those quantities will be identified as the "transverse 
momentum dependent parton distribution functions". They are free from all rapidity divergences, 
including the mixed ultraviolet/non-ultraviolet ones, their evolution is governed by a homogeneous 
differential equation and no auxiliary parameters are needed to resum large logarithms. 



Introduction: Transverse momentum dependent 
(TMD) factorization theorems are essential theoretical 
tools in modern high-energy colliders. High-energy scat- 
tering reactions such as semi-inclusive deep- inelastic scat- 
tering (SIDIS), the (^-dependent spectrum for Drell- 
Yan (DY) heavy lepton pair production or the analogous 
gluon-gluon fusion process through proton-proton colli- 
sions at the LHC, all ask for a TMD factorization theo- 
rems to be established before carrying out phenomenol- 
ogy. Such theorems have mainly dual purpose. First 
is to disentangle long and short distance physics while 
keeping the TMD manifest. This serves to consider qr- 
dcpendcnt spcctrums in which the product(s) of high- 
energy reactions are produced (or tagged) with a non- 
vanishing transverse momentum. Those events consti- 
tute a major bulk of the LHC data. Secondly, when the 
spin is invoked with polarized targets and/or products 
then the nucleoli's three-dimensional structure as well as 
its momentum and spin distributions can be unraveled 
to a maximal extent. 

Establishing such TMD factorization theorems is much 
more complex than the collinear counterparts (where all 
the TMD is integrated over). The complexity arises since, 
in certain kinematical regimes, the relevant modes that 
capture the long-distance hadronic physics, namely soft 
and collinear 1 (e.g., [3]), have the same invariant mass 
but differ only in their relative rapidities. In such in- 



stances the soft and collinear modes can mix under boosts 
(see the discussion in Sec. VII of [4] and references 
therein), thus the scale factorization becomes a subtle 
issue. 

The TMD parton distribution functions (TMDPDFs) 
where introduced by Collins and Soper [5, 6]. Later on Ji, 
Ma and Yuan [7] re-defined the TMDPDF by subtracting 
a complete soft function from the bi-local collinear quark 
field operator. Since then several works have emerged 
[3, 8-14] with seemingly contradictory conclusions. In 
some places it was even concluded that the notion of 
"TMDPDF" simply cannot be defined at the operator 
level contrary to the integrated Feynman PDF. 

Below we argue that only when the soft function is 
partitioned -in rapidity space- in a correct way and 
combined with the collinear matrix elements, one can 
then separate modes with positive and negative rapidi- 
ties, thus obtaining a "two- hemisphere" -like picture. As 
a consequence, we would be able to define TMDPDFs 
(whether polarized or unpolarizcd) which do not suf- 
fer from any bad features 2 . Mainly they will be free 
from all rapidity divergences including mixed ultravi- 
olet (UV) /non- ultraviolet (nUV) ones. Those diver- 
gences pose a serious problem to the renormalizability of 
such physical objects. Our approach is that such mixed 
UV/nUV divergences should not exist, right from the 



2 To best of our knowledge there is not any single treatment so 
1 Depending on the virtuality p 2 of the modes they are either de- far in which the mixed UV/nUV divergences cancel regardless of 



scribed by SCET-II (p 2 ~ ^qcv^ M or SCET-ijt (p 2 ~ q 2 ,) [2]. regulators or prescriptions. 
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start and before rcnormalization, in the proper definition 
of a TMDPDF. Also, and not less importantly, is the abil- 
ity to recover the integrated PDF from the newly defined 
TMDPDF without reference to any particular choice of 
regulators or prescriptions. 

Outline of the Problem: Let us consider the DY 
kinematics in the center-of-mass frame where the mo- 
menta of the two incoming partons initiating the hard 
reaction are p and p. A general vector v is decomposed 
as: = n ■ + n ■ + v± = v +! ^- + v~ 7 -^ + v±, 
where n = (1,0,0,1), n = (1,0,0,-1). Then, p = 
(Q, 0, 0j_) and p = (0, Q, 0j_) where Q is the virtual pho- 
ton mass. The relevant modes that contribute to the 
DY ^-spectrum are collinear [k ~ Q(1,A 2 ,A)), anti- 
collinear (k ~ Q(A 2 ,1,A)) and soft (k ~ Q(X, A, A)), 
where A ~ qr/Q is small. 

As mentioned earlier the soft modes can become n, n- 
collincar under boosts and vice versa. To separate the 
modes one needs to break the "soft-collinear" boost in- 
variance by introducing rapidity cuts which also serve, in 
perturbation theory, to regularize rapidity divergences, 
i.e., when y = iln|fc + /fc~| — > ±oo. 




FIG. 1: Relevant modes for the factorization of DY qr- 
spectrum. Regions A and C are the pure collinear modes 
in the n and n directions, respectively, and B1+B2 are the 
soft modes, rj stand for generic rapidity regulators ( cut- 
offs) necessary to separate the soft and collinear modes. The 
soft function (B1+B2) is split into two pieces along the "sep- 
arating line" k + = k~ , allowing for the separation of the 
whole space into two regions, i.e., two well-defined TMD- 
PDFs: A+Bl (positive rapidities) and C+B2 (negative ra- 
pidities). The thick dashed line is the bisectrix of S(r/ + ,n~) . 



The hadronic tensor for the g-r-dependent DY spec- 
trum can be "factorized" 3 in impact parameter space 
(IPS) as [3] 



M = H(Q 2 



■V 



(1) 



where we show explicitly the rapidity regulator depen- 
dence. S is the relevant soft function and stand 

for pure (anti)collinear contributions. Notice that J^?L 
and S will depend also on other parameters that regulate 
the genuine infra-red (IR) divergencies which exist also 
in full QCD. Pure collinear matrix elements are calcu- 
lated by first integrating over all momentum space and 
then subtracting the "zero-bin" contribution, which is 
the soft limit of the collinear integrands [4] (see also [15]). 
Generally speaking, this should be done on a diagram-by- 
diagram basis. At operator level one can identify those 
soft contaminations with the soft function itself [16-18], 
however this equivalence might get spoiled for certain 
nUV regulators 4 . 

Thus at the level of the factorization theorem itself, one 
should refrain from subtracting the soft function (since 
this would be based on a regulator-dependent arguments) 
and formulate the relevant theorem in terms of pure 
collinear matrix elements and soft functions-whenever 
they are needed-as in Eq. (1). 

At the operator level one defines the collinear matrix 
element as 

where ^ n = W^T^n- W n and T n , collinear and trans- 
verse Wilson lines respectively, guarantee the gauge in- 
variance of the matrix element [19, 20]. In the above 
expression it is implied that the soft contamination has 
not been subtracted yet. Thus we will refer to it as the 
"naive" collinear matrix element. The soft function rele- 
vant for the ^-dependent DY spectrum is given by 

S = <0| Tr \S^ST] (0+, 0~,y ± ) Is^S?} (0) |0) , (3) 



where S n (x) = Pexp 



ig dsn ■ A s (x + sn) and the 
superscript T stand for soft and transverse Wilson lines, 
respectively. The relevant definitions of the collinear, soft 



In Fig. 1 we have denoted by r) generic rapidity reg- 
ulators that separate soft modes from the n, fi-collinear 
ones. Those regulators should disappear when one com- 
bines all the matrix elements within the factorization the- 
orem, since in full QCD there are no rapidity divergen- 



3 The quotation marks indicate that Eq. (1) should be considered 
as an intermediate step towards the final factorization theorem. 

4 Explicit examples can be found in Eq. (64) of [4], and also in [9] 
where the zero-bin contributions vanish beyond tree-level while 
the soft function does not. 



3 



and transverse Wilson lines for DY and DIS kinematics 
can be found in [3]. 

When the soft function is split symmetrically between 
n and n and combined with the pure collincar matrix el- 
ements as in [3], then the TMDPDF defined in this way 
will depend on the two rapidity regulators rf* 1 . The rea- 
son is that the soft function is not symmetric with respect 
to the "separating line" k + = k~ , which separates nega- 
tive and positive rapidities in Fig. 1, but it is symmetric 
with respect to the thick dashed purple line. 

To illustrate the above points let us consider the TMD- 
PDF proposed in [3], 



^old 



(4) 



The A-regulator was used to regularize all nUV diver- 
gences, namely the IR and the rapidity divergences. This 
is why it remains in Eq.(l) when all the matrix elements 
are combined thus generating the genuine IR of QCD, 
which in turn is washed out by confinement. 
Its 0(a s ) result, before renormalization, is 
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5(1 - X n ) - (1 - X n )[l + ln(l - X n )] 



(5) 



where L T = ln(/i 2 6 2 e 27E /4) and V q / q is the Altarelli- 
Parisi splitting kernel of a quark in a quark. The exis- 
tence of two different values for the A-regulator is meant 
to treat in the most general way the supposedly decou- 
pled two collinear sectors. The first line in Eq.(5) con- 
tains mixed UV/nUV term. The second line is what 
would be the matching (or Wilson) coefficient of the 
TMDPDF onto the integrated PDF after an OPE is car- 
ried out, but it also contains unacceptable dependence 
on the nUV regulator. As it is well-known, Wilson coef- 
ficients should be free from any nUV regulators. The last 
line is simply the integrated PDF. It is clear that only 
with the choice A + = A~ = A the mixed UV/nUV di- 
vergence cancels, also the A-dependence from the Wilson 
coefficient. It could be straightforwardly checked that if 
the latter cancelation does not occur, then one cannot 
obtain the integrated PDF from the TMDPDF. Thus we 
could consider the TMDPDF given in Eq. (4) as an ac- 
ceptable candidate for a TMDPDF with the prescription 
of taking the A 's equal. This observation was mentioned 
explicitly in [10]. However in order to avoid such prescrip- 
tion for the regulators a new definition of the TMDPDF 



should be adopted. We will show below that the pro- 
posed new definition will reduce to the one given above 
when the special case: A + = A~ is considered. 

It should be immediately mentioned that if one adopts 
the definition of the TMDPDF given in [3] then the above 
observation about the (non-) cancellation of UV/nUV di- 
vergencies or the A-dependence in the Wilson coefficient 
does not rely on the use of the A-regulator. In other 
words it is a feature of the definition rather than the 
regulators. If one had used off-shellnesscs [14], or quark 
(and gluon) masses, or the regulator as in [9], then 
the mixed UV/nUV and the nUV regulator dependence 
could only be canceled when taking the same values of 
the regularizing parameters for both collinear sectors. 

However for the general case where rj + ^ this is 
not the case, which sheds doubt on the validity of the 
definition given in Eq. (4). Thus one should strive to 
identify the positive rapidity quanta with one collincar 
sector and the negative rapidity quanta with the other 
one. This will be achieved by splitting the soft function 
into two regions, with no overlap between negative and 
positive rapidity regions (unlike the case of just taking 
the square root of the soft function and combine it with 
each of the two pure collinear sectors). We will show 
next that this splitting is a fundamental property of the 
soft function, Eq. (3), which holds to all orders in per- 
turbation theory. In our arguments we will refer to the 
A-regulator, however those arguments are independent 
of the use of any specific regulator. 

Splitting the Soft Function and Definition of 
TMDPDF: In the kinematical region where Q 3> qr 
Aqcd one can perform an OPE of the result in Eq. (1) 
onto the integrated PDFs as 5 

M = H(Q 2 /fi 2 ) C(x n ,x n ;L T , Q 2 /^ 2 ) 

X 4>n(x n ] A-//1 2 ) (/)n(Xn; A+ / fl 2 ) . (6) 

Since the integrated PDFs, 4> n (n), contain just the n,n 
IR collincar divergences, then with the use of the A- 
rcgulator each one of the PDFs will depend on a single 
In ^j- (In ^jj-) to all orders in perturbation theory. Notice 
also that the matching coefficients H and C cannot de- 
pend on any nUV regulator, i.e., they are A-free which 
is also true to all orders in perturbation theory. 

On the other hand we can express the hadronic tensor 
in terms of matrix elements, as in Eq.(l). To separate 
the modes in rapidity we notice that each one of the pure 



C(z„,3^;L T ,Q 2 /V 2 ) = Qn„; L T ) C£{x n \ L T ) 
in the notation of [3]. 



-2D(L T ) 
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collincar matrix elements depend just on one regulator, 
the one which is relevant for their own sector, while the 
soft function depends on both regulators. 

Now, given that all the IR of QCD is generated by 
the two PDFs, as shown in Eq. (6), and that the pure 
collinear matrix elements depend just on their own regu- 
lator, as shown in Eq. (1) combined with the symmetry 
between n and n in the soft function, we immediately 
derive to all orders in perturbation theory the splitting 
of the soft function, which in terms of a generic rapidity 
regulator is 



S( V + , V -) = yS(r]+,v + )S(v-,V-) 



(7) 



In other words, the soft function S(?y + , rj~) can be writ- 
ten as the geometric mean of S(r] + ,7] + ) and S(r)~, r)~) 6 . 
In general, S(rj, r/) would be the soft function where we 
implement the same rapidity regulator for both n and 
n soft Wilson lines. In the case of the A-regulator this 
translates to using the same value of for both Wilson 
lines. For the i/-regulator introduced in [9] it would be to 
use the same v. For any other regulator the implemen- 
tation of Eq. (7) should be straightforward. 

From a pictorial point of view, S'(^ + , rj~) would be the 
area between the lines rj + and rf~ in Fig. 1. Notice that 
for jj + 7^ the soft function is not symmetric with re- 
spect to the separating line, thus its square root contains 
both positive and negative rapidity soft quanta. On the 
other hand, S(r/~ , rj~) would be the area between the two 
rf lines, symmetric with respect to the separating line. 
Then, it is clear that its square root is "half" of the region 
and contains only positive rapidity soft quanta. Thus 



when we combine J n (v~) with y S(rj~ , rj~) we collect all 
the positive rapidity quanta, then defining properly the 
TMDPDF in the n-direction. Same arguments apply for 
the n-direction. 

Since the soft function can indeed be separated into 
positive and negative rapidities along the separating line 
in Fig. 1, we then define the n-TMDPDF as 



(8) 



and for Fn we interchange r\ with 77 + . 

The pure collinear matrix element was calculated in [3] 



6 Eq.(7) can be also derived using equation (10.113) in 



with the A-regulator, 
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which, as explained before, depends just on A , and the 
soft function is 
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thus establishing Eq.(7) at 0(a s ). Combining the 
pure collinear with Si (A - , A - ) that can be extracted 
from Eq.(10), then the newly defined TMDPDF given in 
Eq. (8) is 
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9/9 y2 



-6(1 - x n ) - (1 - x n )[\ + ln(l - x r 



(11) 



From the above result we notice that there are no 
more mixed UV/nUV divergencies as promised, thus it is 
straightforward to rcnormalize the TMDPDF. Also the 
disappearance of the A-dependence from the matching 
coefficient of the TMDPDF onto the integrated PDF (the 
second line in the previous equation) which allows us to 
obtain the latter from the former by simple integration 
in momentum space. 

Finally, the fact the TMDPDF is free from the rapidity 
regulator gives us an evolution equation, 

d 



drj' 



-F n = ; 



(12) 



allowing us to resum large logarithms as it was done in [3] . 
One more remark is in order. The soft function in Eq.(3) 
does not have any reference to the high scale Q 2 , how- 
ever the Q 2 -depcndcncc in Eq.(10) is generated by the 
A-regulator in subtle way, because the parameters that 
regulate the soft Wilson lines, i.e., S ± (see [3]), are re- 
lated to A ± through the large components of the parton 
momenta. Although this observation depends on the par- 
ticular choice of regulators, the fact is that the Q 2 in the 
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TMDPDF is generated by a product of the large compo- 
nents p + and p~ , and not (p + ) 2 or (p~) 2 , thus one con- 
cludes that the TMDPDF is boost invariant, as it should 
be. In other words, the TMDPDF can be extracted from 
a fixed target experiment and used in a collider one, or 
viceversa. 

Factorization Theorem: Given the definition of 
TMDPDF in Eq.(8) we are able to separate positive and 
negative rapidity modes unambiguously, thus the factor- 
ization theorem in Eq.(l), built from ill-defined matrix 
elements, should be re-written as 



rapidity regulators; for the latter one needs an additional 
set of regulators. 
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AI = H(Q 2 /n 2 ) F n (x n ,b; Q 2 , /?) F n {x n , b; Q 2 , /i 2 ) , (13) 

where the two TMDPDFs are the main hadronic ma- 
trix elements. This factorization theorem is valid for 
Q 3> Qt Aqcd, and it can be refactorized in terms of 
integrated PDFs when qx 3> Aqcd- Notice that this step 
needs each TMDPDF to be matched independently onto 
the corresponding PDF, and this is accomplished only if 
one separates positive and negative rapidity quanta as 
proposed above. This one-to-one correspondence allows 
eventually to recover the PDF from the TMDPDF. This 
passage could not be carried out using Eq. (1) unless the 
A's are set equal since there is no clean separation of 
negative and positive rapidities. This requirement can 
be relaxed with our new definition of the TMDPDF. 

Conclusions and Outlook: We have provided a def- 
inition of a TMDPDF which is free from all rapidity di- 
vergencies including the mixed terms that spoil the renor- 
malization of such quantities. This is extremely impor- 
tant for phenomenological applications of TMDs (quarks 
in the case of DY or SIDIS or gluon TMDPDFs for LHC 
physics) since the anomalous dimensions and the Q 2 - 
resummation kernel D [3], and thus the evolution of in- 
dividual TMDs can be unambiguously determined in the 
sense that no un-physical extra parameters are needed. 
Moreover, the resummation of all large logarithms can 
be obtained by "running" between the different relevant 
scales and by solving a homogeneous evolution equation 
as in Eq. (12). Our results can be readily extended to 
define polarized and unpolarized quark and gluon TMD- 
PDFs which could be considered as a generalization of 
the one introduced in [21], as well as TMD fragmentation 
functions. The definition of TMDPDF given in this letter 
reduces to the one proposed in [3] when the two sectors 
are treated identically We finally remark that our treat- 
ment of the TMDPDF applies for massless and massive 
partons as well, where the quark masses combined with 
an artificial gluon mass serve as IR regulators, but not 
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